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Abstract 

A convergence result for a discontinuous Galerkin multiscale method 
for a second order elliptic problem is presented. We consider a heteroge- 
neous and highly varying diffusion coefficient in L°°(n,Rf^) with uni- 
form spectral bounds and without any assumption on scale separation or 
periodicity. The multiscale method uses a corrected basis that is com- 
puted on patches/subdomains. The error, due to truncation of corrected 
basis, decreases exponentially with the size of the patches. Hence, to 
achieve an algebraic convergence rate of the multiscale solution on a uni- 
form mesh with mesh size H to a reference solution, it is sufficient to 
choose the patch sizes as 0(H\ logttf- 1 )!). We also discuss a way to 
further localize the corrected basis to element-wise support leading to a 
slight increase of the dimension of the space. Improved convergence rate 
can be achieved depending on the piecewise regularity of the forcing func- 
tion. Linear convergence in energy norm and quadratic convergence in 
L 2 -norm is obtained independently of the forcing function. A series of 
numerical experiments confirms the theoretical rates of convergence. 

Keywords: multiscale method, discontinuous Galerkin, a priori error esti- 
mate, convergence 
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1 Introduction 

This work considers the numerical solution of second order elliptic problems 
with heterogeneous and highly varying (non-periodic) diffusion coefficient. The 
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heterogeneities and oscillations of the coefficient may appear on several non- 
separated scales. More specifically, let C M. d be a bounded Lipschitz domain 
with polygonal boundary T. The boundary T may be partitioned into some non- 
empty closed subset (the Dirichlet boundary) and its complement IV := 
r\r.D (the, possibly empty, Neumann boundary). We assume that the diffusion 
matrix A e L°° (fi,Rj^) has uniform spectral bounds < a,/3 < oo, defined 
by 

_ . , . , (A(x)v) ■ v (A(x)v) ■ v 

< a := cssmf mf - — ±-!—L < ess sup sup - — =: /3 < oo. (1.1) 

xGfi ueK«'\{0} V ■ V x 6 !i „ G R<i\{o} v ■ V 

Given / € L 2 (f2), we seek the weak solution of the boundary- value problem 

-V • AVu = f in n, 
u = oiiTd, 
v-AVu = Q onTjv, 

i.e., we seek u g i?^(0) := {v G H 1 ^) \ v\r D = in the sense of traces}, such 
that 

a(u,v):= / AWu-Wvdx= / /wda;=:F(w) for all v € Hj}(T£). (1.2) 

Many methods have been developed in recent years to overcome the lack 
of performance of classical finite element methods in cases where A is rough, 
meaning that A has discontinuities and/or high variation; we refer to [H O 
HU El |6] amongst others. Common to all the aforementioned approaches is 
the idea to solve the problems on small subdomains and to use the results to 
construct a better basis for some Galcrkin method or to modify the coarse scale 
operator. However, apart from the onc-dimcnsional setting, the performance of 
those methods correlates strongly with periodicity and scale separation of the 
diffusion coefficient. 

Other approaches [51 1161 [3] perform well without any assumptions of pe- 
riodicity or scale separation in the diffusion coefficient at the price of a high 
computational cost: in [5j [16] the support of the modified basis functions is 
large and in [3] the computation of the basis functions involves the solutions of 
local eigenvalue problems. 

Only recently in |14j . a variational multiscale method has been developed 
that allows for textbook convergence with respect to the mesh size H, \\u — uh\\ h i/q\ < 
Cf t p/ a H with a constant Cf p/ a that depends on / and the global bounds of 
the diffusion coefficient but not its variations. This result is achieved by an 
operator-dependent modification of the classical nodal basis based on the so- 
lution of local problems on vertex patches of diameter 0(H\ log(if _1 )|). The 
method in [TJ] is an extension of the method presented in [T31 HFJ . 

In this work, we present a discontinuous Galerkin (dG) multiscale method 
with similar performance. The method is a slight variation of the method [8], 
in the sense that the boundary conditions for the local problems are now of 
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Dirichlet type. The dG finite element method admits good conservation prop- 
erties of the state variable, and also offers the use of very general meshes due 
to the lack of inter-element continuity requirements, e.g., meshes that contain 
several different types of elements and/or hanging nodes. Both those features 
are crucial in many applications. In the context of multiscale methods the 
discontinuous formulation allows for more flexibility in the construction of the 
basis function e.g., allowing more general boundary conditions [8]. Although 
the error analysis presented in this work is restricted to regular simplicial or 
quadrilatcral/hexahcdral meshes, we stress that all the results appear to be ex- 
tendable for the case of irregular meshes (i.e., meshes containing hanging nodes). 
We refrained from presenting these extensions here for simplicity of the current 
presentation. Under these assumptions, we provide a complete a priori error 
analysis of this method including errors caused by the approximation of basis 
functions. 

In this dG multiscale method and in previous related methods [HI [8] , the 
accuracy is ensured by enlarging the support of basis functions appropriately. 
Hence, supports of basis functions overlap and the communication is no longer 
restricted to neighboring elements but is present also between elements at a 
certain distance. We will prove that resulting overhead is acceptable in the 
sense that it scales only logarithmically in the mesh size. 

In order to retain the dG-typical structure of the stiffness matrix, we discuss 
the possibility of localizing the multiscale basis functions to single elements. In- 
stead of having 0(1) basis functions per element with 0(H\ log(i? _1 )|) support, 
we would then have 0(| log(-ff -1 )|) basis functions per element with element 
support. The element- wise application of a singular value decomposition easily 
prevents ill-conditioning of the element stiffness matrices, while simultaneously 
achieving further compression of the multiscale basis. 

The outline of the paper is as follows. In Section [2] we recall the dG finite 
element method. Section [3] defines our multiscale method, which is then an- 
alyzed in Section HI Section [5] presents numerical experiments confirming the 
theoretical developments. 

Throughout this paper, standard notation on Lebesgue and Sobolev spaces 
is employed. Let < C < oo be any generic constant that neither depends 
on the mesh size nor the diffusion matrix A; a < b abbreviates an inequality 
a < Cb and a rj b abbreviates a < b < a. Also, let the constant Cp/ a depend 
on the minimum and maximum bound bound (a and /3) of the diffusion matrix 
A fTTTj) , 

2 Fine scale discretization 

2.1 Finite element meshes and spaces 

Let T denote a subdivision of D, into (closed) regular simplices or into quadri- 
laterals (for d = 2) or hexahedra (for d = 3), i.e., = UtetT- We assume 
that T is regular in the sense that any two elements arc cither disjoint or share 
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exactly one face or share exactly one edge or share exactly one vertex. 

Let £ denote the set of edges (or faces for d = 3) of T; £ (0) denotes the 
set of interior edges, £(T), £{^o) and £(T N )) refer to the set of edges on the 
boundary of O, on the Dirichlet and on the Neumann boundary, respectively. 
Let T, denote the reference simplex or (hyper)cube and let T P (T) and Q P (T) 
denote the spaces of polynomials of degree less than or equal to p in all or on 
each variable, respectively. Then, we define the set of pieccwise polynomials 

P P {T) := {v : -> R | VT G T, v\ T a F T G Tl p (f)}, 

with TZp G {V p , Q p }, where Ft ■ T — > T, T G T is a family of element maps. Let 
also IIp(T) : L 2 (Vl) — > P p (T) denote the L 2 -projection onto T-piecewise poly- 
nomial functions of order p. In particular, we have (L±o(T)/)|t = l^l" 1 j T f Ax, 
T G T for all / e L 2 (f2). Note that v G P P (T) does not necessarily belong to 
if^fi). The T-piecewise gradient V r w, with (V t w)|t = V(u| T ) for all T G T, 
is well-defined and V r w G (P p _i(T)) d . 

For any interior edge/face e G £(f2) there are two adjacent elements T~ and 
T + with e = <9T~ n dT + . We define f to be the normal vector of e that points 
from T~ to T + . For boundary edges/faces e G £(r) let z/ be the outward unit 
normal vector of fi. 

Define the jump of v G Pk(T) across e G £(f2) by [v] := u|r- — u It+ and 
define [v] := v\ e for e G £ (r). The average of v G P P (T) across e G £ (f2) 
is defined by {«} := (v\ T - + v\ T +)/2 and for boundary edges e G £{T) by 
M := u| e - 

In the remaining part of this work, we consider two different meshes: a coarse 
mesh Th and a fine mesh Th, with respective definitions for the edges/faces £h 
and £h- We denote the Tfr-piecewise gradient by V#i/ := ^7t h v and, respec- 
tively, VhV '■= ^T h v for the 7^-piecewise gradient . We assume that the fine 
mesh Th is the result of one or more refinements of the coarse mesh Th- The 
subscripts h, H refer to the corresponding mesh sizes; in particular, we have 
H G P {T H ) with H\ T = diam(T) =: H T for all T G Th, H e = diame, for all 
e G £h, h G PoC^), with h\x — diam(T) =: hr for all T G 7/j, and /i e = diame 
for all e G £h- Obviously, h < H. 

2.2 Discretization by the symmetric interior penalty method 

We consider the symmetric interior penalty method (SIP) discontinuous Galerkin 
method [3 [HUH]. We seek an approximation in the space Vh '■= Pi(Th)- Given 
some positive penalty parameter a > 0, we define the symmetric bilinear form 
a h : V h x V h -> M by 

a h {u, v) :={AWhU, Vhv) L ^{n) - ^ fa" ' M)i 2 (e) 
+ ({i/ • AVv}, [u]) L 2 (e) - ^-{[u], [v}) L 2 (e) y 
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The jump-seminorm associated with the space Vh, is defined by 



e6£ h (n)U£h(r D ) 



(2.2) 



while the energy norm in Vh is then given by 



•HU:=(||i4 1/2 V h «| 



L 2 (f2) 



+l«l^) 1/2 - 



(2.3) 



If the penalty parameter is chosen sufficiently large, the dG bilinear form (|2.f |) 
is coercive and bounded with respect to the energy norm (|2.3|) . Hence, there 
exists a (unique) dG approximation Uh £ Vft, satisfying 



We assume that (|2.4j) is computationally intractable for practical problems, 
so we shall never seek to solve for Uh directly. Instead, Uh will serve as a reference 
solution to compare our coarse grid multiscale dG approximation with. The 
underlying assumption is that the mesh Th is chosen sufficiently fine so that Uh 
is sufficiently accurate. The aim of this work is to devise and analyse a multiscale 
dG discretization with coarse scale H, in such a way that the accuracy of Uh 
is preserved up to an 0(H) perturbation independent of the variation of the 
coefficient A. 

3 Discontinuous Galerkin multiscale method 

As mentioned above, the choice of the reference mesh Th is not directly related 
to the desired accuracy, but is instead strongly affected by the roughness and 
variation of the coefficient A. The corresponding coarse mesh Th, with mesh 
width function H > h, is assume to be completely independent of A. To en- 
capsulate the fine scale information in the coarse mesh, we shall design coarse 
generalized finite element spaces based on Th- 

3.1 Multiscale decompositions 

We introduce a scale splitting for the space Vh- To this end, let 11^ := Hi(Th) 
and define Vh '■= IlffV/i = P\{Th) and 



Lemma 1 (L 2 -orthogonal multiscale decomposition). The decomposition 



is orthogonal in L 2 (Q). 

Proof. The proof is immediate, as any v € Vh can be decomposed uniquely 
into a coarse finite element function vh '■= IlffW G Vh and a (possibly highly 
oscillatory) remainder v f := (I — Hh)v <G V f , with ||v||| 2 (f2) = H^lli 2 ^) + 



OLh{uh, v) = F(v) for all v G Vh- 



(2.4) 



V f := (f - U H )V h = {v e V h I U H v = 0}. 



v h = v H ®v l , 




L 2 (f2)- 



□ 
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We now orthogonalize the above splitting with respect to the dG scalar 
product ah] we keep the space of fine scale oscillations V f and simply replace 
Vh with the orthogonal complement of V 1 in Vh- Wc define the fine scale 
projection $ : Vh — > V f by 

a,h($v, w) = ah(v, w) for all w G V i . (3-1) 

Using the fine scale projection, we can define the coarse scale approximation 
space by 

VjT := (1 - $)V H . 
Lemma 2 (a/j-orthogonal multiscale decomposition). The decomposition 

v h = v% s e v f , 

is orthogonal with respect to ah, i-e., any function v in Vh can be decomposed 
uniquely into some function G V 1 ^ 8 plus v f G V 1 with |||«|||^ ~ |||w^ s |||^ + 
|||w f |||^. The functions Vj} 3 G Vjj s and v { G V f are the Galerkin projections of 
v G Vh onto the subspaces V^ s and V f , i.e., 

a h (vH S ,w) = a h (v,w) = F(w) for all w G V H IS ', 
ah(v l , w) = ah(v, w) = F(w) for all w G V f . 

The unique Galerkin approximation u^ s G V^ s of it G V solves 

a h (u^ s ,v) = for all v G V^ s . (3.2) 

We shall see in the error analysis (cf. Theorem [8|) that the orthogonality yields 
error estimates for the Galerkin approximation u^ s G V^} s of (|3.2j) that are 
independent of the regularity of the solution and of the diffusion coefficient A. 
However, the space V^ s is not suitable for practical computations as a local 
basis for this space is not easily available. Indeed, given a basis of Vh, e.g., 
the element- wise Lagrange basis functions {Xt.j | T G Th, 3 ' = L ■ • ■ ,r} where 
r = (1 + d) for regular simpliccs or r — 2 d for quadrilaterals/hexahcdra. The 
space V# s may be spanned by the corrected basis functions (1 — B)\tj, T G 
Th, j = 1, ■ ■ ■ ,r. Although Ar,j has local support suppA^j = T, its corrected 
version (1—^)Xtj may have global support in f2, as (|3.1j) is a variational problem 
on the whole domain f2. Fortunately, as we shall prove later, the corrector 
functions 4>T,j decay quickly away from T (cf. previous numerical results in 
[5] and a similar observation for the corresponding conforming version of the 
method 14J). This decay motivates the local approximation of the corrector 
functions, at the expense of introducing small perturbations in the method's 
accuracy. 

3.2 Discontinuous Galerkin multiscale method 

The localized approximations of the corrector functions are supported on ele- 
ment patches in the coarse mesh Th- 
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Figure 1: Example of a one layer patch w^, two layer patch and a three 
layer patch w|., on a quadrilateral mesh. 

Definition 3. For all T € Th, define element patches with size L as 
uj^p :— int(T), 

w£ := int(U{T' G Th \ T' n w^" 1 7^ 0}), jfc = 1, 2, . . . . 

We refer to Figure I3~!?l for an illustration. 

Wc introduce a new discretization parameter I > £ N and define localized 
corrector functions ^ G V f (w^) := {v G V f | v|n\ w £ = 0} by 

a h(4>T,j' w ) = a h(^Tj,w) for all w G V f (w^). (3.3) 

Further, we define the multiscale approximation space 

V™ S ' L = span{A T)j - ^ | T G T H , j = 1, . . . ,r}, 

The dG multiscale method seeks u^ s ' L G V^ S,L such that 

a h (u™' L ,v) = F(v) for all uG V™ S ' L . (3.4) 

Since V^ S ' L C Vfc, this method is a Galerkin method in the Hilbert space Vh 
(with scalar product ah) and. hence, inherits well-posedncss from the reference 
discretization (|2.4[) . 

Moreover, the proposed basis {Xtj — 4>t j \ T ^ Th, j = 1, ■ ■ • , r} is stable 
with respect to the fine scale parameter h, as we shall see in Lemma [7] below. 

3.3 Compressed discontinuous Galerkin multiscale method 

The basis functions in the above multiscale method have enlarged supports 
(element patches) when compared with standard dG methods (elements). We 
can decompose the corrector functions into its element contributions 

4>T,j = X! ^T.jXT', 

T'eTh:T'Cu>%. 
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where XT' is the indicator function of the element T" gTh- 
This motivates the coarse approximation space 

W™ s - L =span({A T j \T € T H , J = 1, . . • , r} 

U {^jXr' |T, T' e T ff , fc^,i = l,.., r}, 

This space offers the advantage of a known basis with element-wise support 
which leads to improved (localized) connectivity in the corresponding stiffness 
matrix. This is at the expense of a slight increase in the dimension of the space 

dim(W™ s ' L ) « L d dim(V]T L ). 

The corresponding localized dG multiscale method seeks w r ^ s ' L € W^ 3 ' L 
such that 

a h (w™ s > L ,v)=F{v) for all v € W™ s L . (3.5) 
Since V^ S ' L C W™' L C Vh, Galcrkin orthogonality yields 

mi ms.L\\\ ^ mi ms,L|M /o r\ 

\\\u h -w H ' \\\h < \\\Uh - u H (3.6) 

i.e., the new localized version (|3.5[) is never worse than the previous multiscale 
approximation in terms of accuracy. However, it may lead to very ill-conditioned 
element stiffness matrices (cf. Lemma [10] which shows that 4>t j£,T' may be very 
small if the distance between T and T" relative to their sizes is large). 

To circumvent ill-conditioning, one may choose a reduced local approxima- 
tion space on the basis of a singular value decomposition. Since the dimension of 
the local approximation space is small (at most proportional to L d ), the cost for 
this additional preprocessing step is comparable with the cost for the solution 
of the local problems for the corrector functions. 

To determine an acceptable level of truncation of the localized basis func- 
tions, we can use the the a posteriori error estimator contribution of the local 
problem from [8], which is an estimation of the local fine scale error. This pro- 
cedure may additionally lead to large reduction of the dimension of the local 
approximation spaces. 



4 Error analysis 

We present an a priori error analysis for the proposed multiscale method (|3.4[) . 
In view of (|3.6[) . this analysis applies immediately to the modified versions 
presented in Section [3. 31 The error analysis will be split into a number of steps. 
First, in Section 14.11 we present some properties of the coarse scale projection 
operator n# . In Section 14. 2\ an error bound for dG multiscale method u^ s 
from p.2[) (Theorem [5J is shown, whereby the corrected basis functions are 
solved globally. Results for the decay of the localized corrected basis function 
(Lemma [10] and Lemma [TT]) are shown, along with an error bound for the dG 
multiscale method u^ s ' from (|3.4[) (Theorem IT^j) , where the corrected basis 
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functions are solved locally on element patches. Finally, in Section POl we show 
an error bound given a quantity of interest (Theorem I14p , leading to an error 
bound in L 2 -norm ( Corollary [T5)) . 

We shall make use of the following (semi)norms. The jump-seminorm and 
energy norms, associated with the coarse space Vh, are defined by 

\*\%-= E ]^IIMIl!»(e). 

eee H (n)u£ H (r D ) e 

|||. Ilk ^(ll^ 1/a Vir.|li.( n ) + l-l!r) 1/a , 

respectively, along with a localized version of the local jump and energy norms 
(|2.2|) and ()2.3j) on a patch u> C f2, where u is aligned with the mesh Th, given 

by 

\ E fllMlli»(e)» 

l||.|IU, w ^(P^.n^M + l-lL) 172 - 

The shape-regularity assumptions Ht ~ ^ e for all e 6 <9T : T e Th and Ht ~ -He 
for all T e <9T : T e Th will also be used. 

4.1 Properties of the coarse scale projection operator 11^ 

The following lemma gives stability and approximation properties of the oper- 
ator Hh- 

Lemma 4. For any v € Vh, the estimate 

H~ l \\v - U H v\\ L 2 {T ) < o^^HHII^T, 

is satisfied for all T G Th ■ Moreover, it holds 

p- 1/2 \\\TlHv\\\H + \\H-\v - njrt;)|| ia( n) < « _ 1/2 1 1 1 ^1 1 U ■ 

Proof. Theorem 2.2 in [T2], implies that for each v £ Vh, there exists a bounded 
linear operator XJj : Vh — > Vh fl H 1 ^) such that 

^ii^Vh^-^h^^ + II/i- 1 ^-^)!!^^) <a- 1/2 H h . (4.1) 

Split v = v c + v d G Vh into a conforming, w c = T^y, and non-conforming, 
u d = v — part. We obtain 

h-'Wv - n H v\\ L 2 {T) < H-Wiv* - n H v c \\ L 2 (T) + \\v d - rW|| L2(T) ) 

< ||V h «|| + ||V fc (« - v c )\\ L 2 (T) + H-^v^t)) (4.2) 

< a _1/2 |||w||U,T 
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using the triangle inequality, stability of the L 2 -projection, and (|4.ip . Also, 

^\\[v c -IL H v}\\l Ke) 



\n H v\\\ 2 H = Y \\VAv(u H v-u (T H )v)\\ 2 L2iT) + Y,^\\[yc-W 2 



teTh eer 
% E P[^\\ v -MT H )v\\h (T) + ^\\v c -IL H v\\l HT) 

T£Th ^ 

<r 2 infill 2 

using the triangle inequality, (|4.1[) , and (|4.2|) which concludes the proof. □ 

The operator Tin is surjective. The next lemma shows that, given some vh G 
Vh in the image of Hh there exists a H 1 -conforming pre-image v € HJ^vh C Vh 
with comparable support. 

Lemma 5. For each Vh G Vh, there exists a v G Vh H i? (0) smc/i i/iai Ilfl-u = 
IIMIU ^ C' (9 / Q |||ui/|||i/, and supp(v) C supp(v H ). 

Proof. Using (|4.1[) but on space Vh gives, for each v G V#, there exists a 
bounded linear operator 1^ : Vh — > Vi? n H 1 ^) such that 

/T 1/2 P 1/2 V ff ( V - Z&tOllwp,) + ll^" 1 ^ - 2&«)||i»(n) < a- 1/2 \v\H- (4.3) 
We define 

v := 2 c H v H + Y ( v h( x j) - x h v h{x 3 )) 9 T ,j, 

TeTh.. j=l,-,r 

where 9tj 6 Vj, fl Hq(T) arc bubble functions, supported on each element 
T, with U H G T ,j = A T j and |||0 T ,j||U < j3(T)H d ~ 2 . Observe that supp(w) C 
supp(w^r). The interpolation property follows from 

U H v = l c H v H +U H 2J ( v h(xj) ~ ZHV H {xj)) 6j, 

TeT H . j=l,...,r 

= T h v h + v H - T-hVh = v H - 
To prove stability, we estimate |||w|||/t as follows: 

\\\v\\\ 2 <\\A 1 / 2 Vlj 1 VH\\h {n) + Y \vh(x 3 )-2 c h vh(x 3 )\ 2 \\\9 3 \\\ 2 

T£Th, 3=1,. ..,r 

< WA^hThVhW 2 ^ +P\\H- 1 {v h -ZHVH)\\h {U) 
<C 2 /a \\\v H \\\ 2 H 



using the inverse estimate ||v||l<=°(t) < H d / 2 |M|,L 2 (T) f° r all v € Vh, and the 
estimate 114.31). □ 



Remark 6. Note that 6» TJ G V h ni/o( T ) for all T e Tff (fulfilling the conditions 
in Lemma [5j can be constructed using two (or more) refinements of the coarse 
scale parameter H . We can let 9t.j G Vh' CiHq(T) where Vh' C Vh and h < h' < 
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2~ 2 H. This does not put a big restriction on h since the mesh Th is assumed 
to be sufficiently fine to resolve the variation in the coefficient A, while the 
parameter H does not need to resolve A. 

The following lemma says that the corrected basis is stable with respect to 
the fine scale parameter h in the energy norm (|2.3|) . this is not a trivial result 
since the basis function {\tj\T G Th, j = 1, • • • , f} are discontinuous. 

Lemma 7 (Stability of the corrected basis functions). For all, T G Th, j = 
1, . . . , r and L > G N, the estimate 

II I At, j - ^TjllU ~ C/3/alll'Vrj|||ff, 
is satisfied, independently of the fine scale parameter h. 

Proof. For any T G Th, j — lj • • • > f, by Lemma [5] there exists a b such that 
v = X T ,j - be V[{uj^), and |||6||| h < Cp /a \\\\ T ,j\\\H- We have 

1 1| At j - 0Tj|llh ~ a h (X T j - 0tj, A T j - 4>tj) = a h(^T,j - <f>T,j: ^tj - v), 

< 0/i(Atj - </>Tj; ~ ~ ^T,j II Ul II ^T,j 1 1 Iff, 

which concludes the proof. □ 
4.2 A priori estimates 

The following theorem gives an error bound for the idealized dG multiscale 
method, whereby the correctors for the basis are solved globally. 

Theorem 8. Let u h G V h solve and u^ s G V% s solve (03), then the 

estimate 



,ms I 



|||u* - UiTlllfc < Cia-^WHif - n H f)\\ L 2 {nh 
is satisfied, where C\ neither depends on the mesh (h or H) size nor the diffusion 
matrix A. 

Proof. Let e := Uh — u 7 ]} s =Uf G V f , then 

IINIIh ^ a h (e,e) = (f,e) L 2 {n) = (/ - U H f,e - H H e) L 2 {n) 
< \\H(f - nHf^^WH^ie - H H e)\\ 2 L2{n) 

<4=||^(/-n H /)|| L3( Dll|e||| h> 

using Lemma [2] Lemma [T] Cauchy-Schwarz inequality, and Lemma |4l respec- 
tively. □ 

Definition 9. The cut off functions Cfp D G PoiTh) are defined by the conditions 

,d,D } _ 1 
St — 

Ct' D |o\w^ = 0, 
\\[<WU-mSjp^ forallTGT,, 
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and that £ T ' is constant on the boundary d(uj^ \ w^). 

The next lemma shows the exponential decay in the corrected basis, which 
is a key result in the analysis. 

Lemma 10. For all T € Th , J ' = 1, ■ • • , f, the estimate 

lll(A, - <h ti ) - (Aj - <fc d )\\\ h = \\\<hj ~ <&J\\h < Csl L \\\<pT,j - A||U, 

is satisfied, with C 3 = CCp,, < 7 < 1 given by 7 := {p2.)-^r , C 2 = C'C^ a , 
and L = kl, k, I > 2 € N, noting that C and C are positive constants that 
are independent of the mesh (h or H ), of the patch size L, and of the diffusion 
matrix A. 

Proof. Define e := — 4>t j = 4>T,j — ¥t y We nave 

lll e lll^ iS ah(e,<h,j - <t>T,j) = a-h{e,<h,j ~ v) < \\\e\\\ h ■ \\\(f>T,j - v\\\ h , (4.4) 
for v £ V^(w^ fc ). Let ( := (!^~ 1 ' ek ; then by Lemma [5] there exists a & such that 

v = {fa - b e v£K fc ), n ff 6 = iLjC^tj, |||&||U ;$ qs/allinjrC&ylllir. and 

supp(6) C supp(II//C0Tj")- Then, we have 



< 1 1 1 0r i3 - - C0t,j \\\h + Cp /a \\\n. H {C<h,j - <frj )\\\h 

<cf 3/a \\\^~CM\\ h . 



(4.5) 



Furthermore, using the properties of £ we have 

llVIVfeCl-Cfell^cfi),^ \\^V h <j>T,j\\ L *{n\^y ( 4 - 6 ) 

and 



l(i-cfci^= £ f ll[(i-cfc]|i! 2(e) 

eG£(0)U£(r D ) 6 

^ E f (l^ 1 - C}[^]||| 2(e) + ||fe}[l - C]||| 2 (e)) 



ee£(n)uf(r D ) 

< E (^lfe]|l! 2(e) + £^ll{^}|li 2(e ) 



ee£(n)U£(rr,): 



(4.7) 



E flfe]||!2 (e) + ^II^tj -n^ TJ || 



2 



< 



enn\u e T k - Vo 

^2 



c '^/alll^r,illU,n\a>"- 1 ' 
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using a trace inequality and Lemma 01 respectively. Combining (|4.4[) . (|4.5[) , 
(f4~6|) . and (g3) yields 



C^IIU^^/JII^IH,^-,. (4.8) 



To simplify notation, let m := £(k — 1) — 1 and M := £k — l. For r\T := 1 — CJt 
we obtain 



IfeHI^M < llhr^HI^ < a h (jiT<hj,m"i>T,i), (4-9) 

where 

+ ( _2 ({^ • A VriT<t>T,j}, [riT<f>T,j]) + ^([?7T0Tj], [??T</>Tj]) 

eG£(n)uf(r I3 ) ^ e 

(4.10) 

For the first term on the right-hand side of (|4.10p . we have 

(AVhVT<t)Tj,VhVT(f>T,j)L2(n) = (AVh<t>T,j, ^hVT^T,j)L^(Q), (4-H) 

since T]t is constant on each element T e 7ft,; for the other terms we use (|A.3|) 
and (|A.4[) (with v — r/x, w = v ■ AXi '4>T.j and u = 4>r.j)- We can, thus, arrive to 

Ill^jlllftUv^ - a h(VT<pT,j,VT(f>T,j) = Cb h {(j)T,j,riT<h,j) 

+ Y, (V2({^ • AVcj> T ,j}, [VT} 2 [<t>T,j})me) - l/4([^ • AVfoj], {VT} 2 {^ 3 })LHe) 

- ^(M 2 , [<M V(e) + ^(M 2 , {0T J } 2 )L 2 (e)), 

(4.12) 

using gU), (|4~T0|) . and (j4JT|) . Note that, 

£ (l/2({l/ • AV^i}, W'^-Dl^) - l/4([i/ • AV^Tj], M 2 {^T,i})i=(e) 

- ^([^] 2 , [<MV( e) + ^(M 2 ,fe} 2 )^(e)) 



< 



eef(O): 



+ ^(ll[0T. J ]|li 2(e) + !|{0T J }|li 2(e) )) 



~ (^J^2"H j4V '? !, Tj!|l2(T+UT-)!I'/ , TjI|l^T+UT-) + ^H^H^'Hi^T+UT-)) 

ee£(0): T T 

</?r 2 !|F T 1 (0 Tj -n^ TJ )||2 2(<Va)??) < ^ /Q r 2 |||0 Tj |||L-\^- 

(4.13) 
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Now we bound the term 

ah{(t>T,j-,'rjT ( t>Tj) = dh((t>T,j,VT^T,j ~ b ) + a h{<t>T,j , &) = 0>h((f>T,j,b) 

^ III^T,j||U,a,M\ w?? |||6|||^ w M Va;?; (4.14) 

< C/3/a 1 1 1 1 1 U.uM^m 1 1 1 IIjj 7? T </) T , j 1 1 1 fl>M\ w m . 

Furthermore, we have that 



= \\VA\7 H n H {rt - U (T h )Vt)<I>t,j\ 



2 



(4.15) 



eg g(n)ug(r D ): e 

< ^hh- 1 ^^ - n (rH)r?|)^iii 2(W M VS?) 

< pi-*\\H-\4> Td - n^ T ,^||£ 2(<w?) 
<^- 2 |||^|||^^, 

using a trace inequality, inverse inequality, and Lemma 21 respectively. Com- 
bining the inequalities (|4.12j) , (|4.13|) . (|4.14l) . and f|4.15[) yields 

HlfcylllU< ^ ylll^lllWv? ^ Ylll«rjlllU«?- 

where C-2 = C'C^ a . Substituting back to £ and k and using a cut off function 
with a slightly different argument, yields 

<-<(t)*" 1 iii^iii^\^- 

which concludes the proof together with (|4.8|) . □ 
Lemma 11. For a/Z, T G 7jj, j = 1, . . . , r, £/ie estimate 

HI £ ^fe-^OIII^c^ £ KHll^-^-lll 2 , 

TeTh, j=l,—,r TeT H , j=h-,r 

is satisfied, with C4 = CC| / Q and C positive constant independent of the mesh 
(h or H), of the patch size L, and of the diffusion matrix A. 

Proof. Let w = Y^teTh j=i r v j( c t ) T.j — 4>t j)i an d note that 

ah((/)T. 3 - Ar,j, w - Ctw + b) = 0, 

r (4.16) 
ah(4>T 3 ~ *>T,j,w- Ctw + b) = 0, 
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where Ct := C^ +1 ' i+2 , using Lemma [S] and the property of the cut-off function. 
From ((¥7T6|) follows that a h (w, (1 - ( T )w + b T ) = for all w E V i . We obtain 

III E v 3 ( y ( t>T,j - 4>T,j)\\\h< E Vja h ((t)T,j -<t>r,j,w) 
T€Th, j=l r T&Th, 3=1 r 

= E V 3 a h(<pTj ~ (j>T,j,CTW - b) 
TeT H . j=l,-,r 

< E l^l-lll^-^-IIUdllCrHIU + IIHIW (4.17) 

T£Th. j=l,-,r 

< E l^l-lll^-4 J IIU(IIICT«;||U + ^ /a |||n i? CTHIk) 

T£Th, j=l,—,r 

< E l«il-IH^-^rjlllhC?/alllCTHIIfc 

TeT H , j=l,-,r 

From (|4.6p and (|4.7p . we have 

IIICHIU = IHCrwIIU ,^+2 < C/3/al||^|IU,4+=- ( 4 - 18 ) 

Then, further estimation of f)4. 1 7[) can be achieved using (|4.18j) and the discrete 
Cauchy-Schwarz inequality: 

ill E v A<hj - <&,j)\\\h 

T<ETh, j=l,-,r 

/ \ 1/2 / XV2 

<c%A e N 2 ife-0yd E m-iiiU- 

\TgTh, j=l,...,r J \TeT H , 3=l,-,r 

<cl /a L^. ( E KflH^-^-lll^ ■IIHIU- 

\TGTh, j=h-,r J 

Dividing by w on both sides concludes the proof. □ 

The following theorem gives an error bound for the dG multiscale method. 

Theorem 12. Letu g H^(Cl) solve $H2§ andu™ s ' L G V^ s L solve (03). Then, 
the estimate 

WW - u™' L \\\ h <\\\u - u h \\\ h + da-^WHif - U H f)\\ L2{n) 
+ C 5 \\H- 1 \\ L ^ (n) L d l 2 1 L \\f\\ L , m , 

is satisfied, with < 7 < 1, L from Lemma \1(A C\ from Theorem^ C5 = 
^^S/a^t^ 2 ^ w ^ ere C3 ! ' s from Lemma \10\ and C4 is from Lemma C a 
positive constant independent of the mesh (h or H), of the patch size L, and of 
the diffusion matrix A. 
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Remark 13. To counteract the factor 1 1 1 1 1 =° ( ) i n the error bound in Theo- 
rem[l2l we can choose the localization parameter as L = [Clog(| — 1 1 1^=° (o) )1 • 
On adaptively refined meshes it is recommended to choose L = [Clog(if -1 )] . 

Proof. We define u^ s ' L := J^teTh j=i r u H^T( x j) ( l ) T,j- Then, we obtain 

Ml ms,L|M ^ in -ms.Liii 

|| |w- u H \\\h < \\\u - U H \\\h 

^ 1 1 1 III i 1 1 1 tyis ill i ill ms ~ms,L, \ \ \ 

< \\\u - Uh\\\h + \\\Uh - u H \\\ h + \\\u H -u H \\\ h 

< \\\u - Uh\\\h + \\\Uh - Ux S \\\h + HI E u H^T( x j)( ( t ) T,3 - (l>T,j)\\\h- 



(4.19) 



Now, estimating the terms in (|4. 19[) . we have 

UK - «£ S IIU < c ia -^ 2 \\H(f - n ff /)|| £ . ( n), 

using Theorem [51 and 

in e u^ixj^j-^M 

TeTh, j=1,-,r 

<C A L d Y, \ u T: T {^)\ 2 Mt, 3 - (4.20) 

TeT H ,j=l,-,r 

< c 4 c!L d ^ L E I^M^Hfe - \j\wl 

TeTh, 3=1,. .-,r 

using Lemma lTTl and Lemma fTOl respectively. Further estimation, using LemmaO 
yields 

E K.V(^)| 2 IH^"A TJ |||2 

TeTh , j=l,...,T 

^ C l/ a E I«S.t(^)I 2 II|At,,|||^ 

T&Th, j=h-,r 

~ c l/ a P E \ u h s ,t( x 3)\ 2h t 2 \\^tA\ 2 l*(t) (4.21) 

TeTh, j=l,-,r 

= C l/aP E \\ H T 1u ™t{. x 3)^Tj\\ 2 L z(T) 
TeTh, j=l,-,r 

<Cj,J\\ E HT^A^Tjhw- 

TETh , 3=1,. ..,r 

Furthermore, using a Poincare-Friedrich inequality for piecewise H 1 functions, 
we deduce 

II E H T lu H*T( x j)^Tj\\L2(Q) 
T£Th, 3=l,...,r 

<|| E H^u^ T { Xj )n H {X T>j - fhMhw U22) 

TET H ,j=l,...,r 1 ' ' 

<oT x \\\U- x u^\\t 
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□ 



Combining (|i^DJ) . (|4~2"Tj) and 022), we arrive to 

lll«r-«r ,L |ll^^/a^ 1/a C'3||^- 1 |Uoc (n) L«'/Vll/IU»(n). 

4.3 Error in a quantity of interest 

In engineering applications, we are often interested in a quantity of interest, 
usually a functional g £ L 2 (tt) of the solution. To this end, consider the dual 
reference solution (|2.4[) : find </>h £ Vt such that 



a/i(u, 4>h) = g(v) for all v £ V fc , (4.23) 
and the dual multiscale solution (|3.4[) : find 0^ h £ h such that 

ah(v,<t%% L )=g(v) for all « G V£ >fc . (4.24) 

Theorem 14. Lef u £ sofoe (fO]l . w" s ' i £ V# 8 ' L so/we ([g^]) . and let 

g £ L 2 (f2) be the quantity of interest. Then, the estimate 

\g(u) - g{u^ h )\ < \g{u) - g{u h )\ + \\\u h - i% s,i |||fe|||(?f>/» - <j>™ ,L \\\h, 
is satisfied. 

Proof. From (|4.23l) and (|4.24[) . we obtain the Galcrkin orthogonality 

a h {v^ h -4>T L ) =0 for alH- £ V™ S ' L . (4.25) 
Using the triangle inequality, we have 

\g(u) - g{u™ s ' L )\ < \g(u) - g(u h )\ + \g(u h ) - g{u™ s ' L )\- 
Finally, observing that 

\g(u h - ul H )\ = \a h (u h - u^ H ,<f> h )\ 

= \a h (u h - Uh iH ,<f>h - <f>H,h)\ 

£ I IK -^.hIIUIIK - ^hmWU, 

using (|4.25|) . concludes the proof. □ 

Corollary 15. For g{v) = (it/, — ujj h , v)i,2/q-\, the following L 2 -norm error 
estimates hold: 



and 



I L II ^ II ll i ill ms,L 1 1 1 1/2 1 1 1 , ,ms,L|||l/2 

\u-u Hth \\ L 2 {n} < \\u -u h \\ L 2 {n) + \\\u h -u H HI,/ \\\<ph-4>H \\\h » 

\W - "H,fclU 2 (n) S \\ u - uh\\mn) + H\\\u h - u™' L \\\h, (4.26) 

for L = [Clog(i? -1 )] with C sufficiently large positive constant independent of 
the mesh parameters. 
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Remark 16. As expected, if we are interested in a smoother functional, a 
higher convergence rate is obtained for \g(v,h — u ff S ' L )|- For example, given 
the forcing function for the primal problem / £ if m (Tjf), a quantity of interest 
g € H u (Th) (with H°(Th) denoting the standard L 2 (VL) space), and choosing 
L = |~Clog(-ff -1 )] with large enough C, gives 

\g(u-uT L )\ < \g{u)-g{u h )\ + H 2+m+n { \f\ H ^ T) )( £ \g\ HH T))- 

TETh TeTh 

5 Numerical Experiments 

Let f2 where be an L-shapcd domain (constructed by removing the lower right 
quadrant in the unit square) and let the forcing function be / = l+cos(27ra;) cos(27ry) 
for (x,y) £ £1- The boundary T is divided into the Neumann boundary IV := 
T n ({(x 7 y) : y = 0} U {(x, y) : x = 1}) and the Dirichlet boundary = 
r \ r^v. We shall consider three different permeabilities: constant A\ = 1, 
A-2 = A2(x), which is piecewise constant with periodic values of 1 and 0.01 with 
respect to a Cartesian grid of width 2~ 6 in the x-direction, and A 3 = A 3 (x 7 y) 
which piecewise constant with respect to a Cartesian grid of width 2~ 6 both 
in the x- and y-directions and has a maximum ratio (3 /a = 4 • 10 6 . The 
data for A3 are taken from layer 64 in the SPE benchmark problem, see 
http://www.spe.org/web/csp/. The permeabilities A2 and A3 are illustrated 
in Figure [2] For the periodic problem many of the corrected basis functions 




(a) (3/a = 10 2 (b) /3/a « 4 ■ 10 6 



Figure 2: The permeability structure of A2 (a) and A3 (b) in log scale. 

will be identical. For instance, all the local corrected basis in the interior are 
solved on identical patches, thereby reducing the computational effort consid- 
erably. In the extreme case of a problem with periodic coefficients on a unit 
hypercube, with period boundary conditions, the correctors 4>T.j, j = 1, . . . ,r, 
will be identical for all T E Th- 

Consider the uniform (coarse) quadrilateral mesh Th with size H = 2~ l , 
i = 1, . . . , 6. The convergence rate — p/2 corresponds to 0{H P ) since the number 
degrees of freedom sa H~ 2 . The corrector functions p.3[) are solved on a subgrid 
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of a (fine) quadrilateral mesh Th with mesh size 2 -8 . The mesh Th will also act 
as a reference grid for which we shall compute a reference solution Uh, £ Vh (|2.4p 
on. Note that the mesh Th is chosen so that it resolves the fine scale features 
for A4, i = 1, 2, 3, hence we assume that the solution Uh is sufficiently accurate. 

5.1 Localization parameter 

If / G H m (T H ) we have the bound 

\\H{f-HHf))\\mn)< E H k+1 \f\ Hk(T)l (5.1) 
TeT H 

where k = for m = 0, k = 1 for m = 1, and k = 2 for m > 1. Hence, to 
balance the error in between the terms on the right-hand side of the estimate in 
Theorem[12l different constant C has to be used for the localization parameter, 
L = \C\og(H~ 1 )~\, depending on the element- wise regularity of the forcing 
function / on Th- Figure [3] shows the relative error in the energy norm 




N dof (#Degrees of freedom) 

Figure 3: Diffusion coefficient A\ = 1. Relative energy-norm error against -/Vdof, 
for different values of C for the localisation parameter L. 

\\\ u h — ul H S,L \\\h/\\\uh\\\h and FigureHlthe relative error in the L 2 -norm \\uh — 
u^ s ' L \\L 2 (fi)/\\ u h\\L 2 {n) between Uh and u^ s ' L against the number of degrees of 
freedom N dof « 0(i/" 2 ), using different constants C = 1,3/2, 2, 5/2. With the 
chose C = 5/2, the errors due to the localization can be neglected compared to 
the errors from the forcing function, both for the energy- and for the L 2 -norm. 
For / i? 1 (7h), C — 3/2 is sufficient since (|5.1[) gives linear convergence. In the 
following numerical experiments we shall use C — 2, since this value seems to 
balance the error sufficiently. Note that the numerical overhead increases with 
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10 1 10 Z 10 3 10 4 10 5 

N (#Degrees of freedom) 



Figure 4: Diffusion coefficient Ai = 1. Relative L 2 -norm error against N^of, for 
different values of C for the localisation parameter L. 

C as the sizes of the patches T e 7~h, increases with L — \Clog(H~ 1 )~\ . 
This results in both increased computational effort to compute the corrector 
functions and reduced sparseness in the coarse scale stiffness matrix. 

5.2 Energy-norm convergence 

Let the localization parameter be 

L = [~2 log(i7 -1 )] . Figure [5] shows the relative error in the energy norm plotted 
against the number of degrees of freedom. The different permeabilities Ai, 
i = 1,2,3, and the singularity arising from the L-shaped domain do not appear 
to have a substantial impact on the convergence rate, which is about 3/2, as 
expected. We note in passing that using standard dG on the coarse mesh only 
admits poor convergence behaviour for A^ and for A3. This is to be expected, 
since standard dG on the coarse mesh does not resolve the fine scale features. 

5.3 L 2 -norm convergence 

Again, set L = \2 log(£f -1 )] . Figure|B]and Figure[?l shows the relative L 2 -norm 
error again the number of degrees of freedom between, Uh and h and between 

u h andnf/uf^, viz., -n_ffU^ J5 ' L || L 2( f2) /||u /l || i 2( a) , respectively. In Figure [6] 

we see that the i 2 -norm error between Uh and uJ^ s ' L converges at a faster rate 
than in the energy norm (convergence rate —2 compared to —3/2, respectively,) 
as expected from (|4.26[) . In Figure [7] only the coarse part of u^ s ' L is used (i.e. 
nff% s,i ); nevertheless it appears to have a faster convergence rate than —1/2, 
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10 1 10 Z 10 3 10 4 10 5 

N (#Degrees of freedom) 



Figure 5: Relative energy- norm error against iVdof, for C = 2 in the localisation 
parameter L for the the diffusion coefficients A\ , A2, and A3 . 

except for the case of the permeability A3. 




10 1 io z 10 3 10 4 10 5 

N (#Degrees of freedom) 



Figure 6: Relative i 2 -norm error against TVdof, for C = 2 in the localisation 
parameter L for the the diffusion coefficients Ai, A2, and A3. 



22 




Figure 7: Relative L 2 -norm error against iVdof, for C — 2 in the localisation 
parameter L for the the diffusion coefficients A±, As, and A3. 
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A Equalities for averages and jump operators 

We derive equalities for averages and jump operators across interfaces where 
the functions v and w have discontinuities. Using [vw] — {v}[w] + [v]{w} and 
{«}{«;} = {vw} — l/4:[v][w], we have 



{ot}[wi[ ={ui}{u}[uu] + l/4[v][u>][vu] 

={w}[v 2 u] — {w}[w]{w?i} + l/4[w] [w] [vu] 
~{w}[vu] — [v]{u;}{u}{ii} — l/4[v]{w}[w] [u] 
+ l/A[v] 2 [w]{u} + l/4[v]{v}[w][u] 



(A.1) 
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(A.2) 



and 

{OTu}[zm] ={v}{iw}[it] + 

={v 2 w}[u] — 1 / 4[v}[vw][u} + {vw}[v]{u} 
={v 2 w}[u] - l/A[v] 2 {w}[u] - 1/A[v]{v}[w][u] 
+ {v]{v}{w}{u} + l/4[v] 2 [w]{u) 
Combining (|A.1|I and ()A.2j) wc obtain 

2{vw}[vu] = {w}[v 2 u] + {v 2 w}[u] + l/2[v] 2 [w]{u} - l/2[v] 2 {w}[u] (A.3) 

Also, 

[tm] =[u]{t)}[t)u] + [«]{u}[!jk] 

= [u][w 2 u] — + [jj]{m}[wu] 

= M[w 2 u]) - [«] [«]{»}{«} - 1/4MMHM (A.4) 

+ [v}{u}[v]{u} + [v]{u}{v}[u] 
= [u][v 2 u] - l/4[v] 2 [u} 2 + [v] 2 {u} 2 



